Structure and dynamics of interfaces between two coexisting liquid crystalline phases 
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The phase-field-crystal model is used to access the structure and thermodynamics of interfaces 
between two coexisting liquid crystalline phases in two spatial dimensions. Depending on the model 
parameters there is a variety of possible coexistences between two liquid crystalline phases including 
a plastic triangular crystal (PTC). Here, we calculate numerically the profiles for the mean density 
and the nematic order tensor across the interface for isotropic-PTC and columnar-PTC respectively 
smectic A-PTC coexistence. As a general finding, the width of the interface with respect to the 
nematic order parameter characterizing the orientational order is larger than the width of the mean 
density interface. In approaching the interface from the PTC side, at first the mean density goes 
down and then the nematic order parameter follows. The relative shift of the two profiles can be 
larger than a full lattice constant of the plastic crystal. Finally, we also present numerical results for 
the dynamic relaxation of an initial order-parameter profile towards its equilibrium interfacial profile. 
Our predictions for the interfacial profiles can in principle be verified in real-space experiments of 
colloidal dispersions. 

PACS numbers: 68.08.De, 61.30.Dk, 64.70.D-, 82.70.Dd 



I. INTRODUCTION 

Liquid crystals typically composed of anisomet- 
ric molecules or colloidal particles form interesting 
mesophases, which are neither completely liquid nor crys- 
talline yy. The simultaneous presence of translational 
and rotational degrees of freedom gives rise to phases, 
which exhibit a different degree of ordering for the trans- 
lational and orientational order [2l, W . Rotator solids or 
plastic crystals, for instance, are translationally ordered 
but orientationally disordered, while nematics, on the 
other hand, possess orientational order in the absence 
of translational order. Clearly, there is also the fully 
disordered isotropic phase and the fully ordered crys- 
talline phase, but there are even more intermediate liquid 
crystalline phases (such as, for example, smectic- A and 
columnar phases) with different degrees of translational 
and orientational order that are stable for appropriate 
thermodynamic conditions. 

It is a formidable task of statistical physics to predict 
the existence and stability of the different liquid crys- 
talline phases for a given interaction as a function of mean 
density and temperature. This has, for example, been 
done by computer simulations of simple model systems 
[4, 5 and also by molecular density functional theory 
[6, JJ and more phenomenological approaches [8^ . Typ- 
ically, phase diagrams of liquid crystals exhibit regions, 
where two phases of different kind of ordering coexist. At 
equal pressure, chemical potential, and temperature, co- 
existence implies that there is a stable interface between 
the two coexisting phases. In mean- field theories (which 
neglect interfacial capillary wave undulations), the inter- 
face has a characteristic width of typically several parti- 
cle sizes and exhibits profiles of the mean density and of 



the degree of orientational order depending on the spa- 
tial coordinate perpendicular to the interface. For the 
liquid-solid interface, see, for example, Refs. P4T4]. 

While there has been a large effort to explore the 
gas-liquid [15 and liquid-solid [16 interfaces of spher- 
ical particles (see also Refs. [17-23 ), much less effort 
has been devoted to the particle-resolved structure and 
thermodynamics of the interface between two coexisting 
liquid crystalline phases. Extensive studies have been 
performed for the isotropic- nematic interface, which has 
been accessed by experiment, computer simulation, and 
theory (see, for example, Refs. ^24-32 ), but there are 
fewer considerations of the isotropic-smectic [33H36] and 
the nematic-smectic interface [37]. However, to the best 
of our knowledge there is no investigation of an inter- 
face, where one of the coexisting phases is plastic or fully 
crystalline. This is, of course, a nontrivial task, since 
there is a complex dependence of the interface structure 
on the (relative) orientation of the two phases. Even for 
the isotropic-crystal coexistence there is a complex orien- 
tational dependence culminating in Wulff 's construction 
for the equilibrium crystal shape [23^. Nevertheless, it is 
important to have information about the interface, since 
nucleation and growth phenomena of a metastable phase 
in a stable phase occur via interfaces [36l |38l |39] . 

In this paper, we close this gap and study liquid crys- 
talline interfaces also for crystalline phases. We use 
a phase-field-crystal (PFC) model, which is a minimal 
model to describe freezing for isotropic particles on the 
molecular (i.e., interparticle) scale [4QH42] and can be 
justified from microscopic density functional theory j42l - 
[44]. The traditional PFC model [40 was later general- 
ized to anisotropic particles in two [45 and three [46] 
spatial dimensions allowing for liquid crystalline phases. 
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The generaUzed theory is formulated in terms of three 
order-parameter fields, namely the reduced translational 
density '^(r), the local nematic order parameter 5'(r), 
and the mean orient ational direction n(r) that is also 
called the "nematic director" . While the traditional PFC 
model [40] has two free parameters, the liquid crystalline 
PFC model in two dimensions [45 has five independent 
couplings. This widely opens the parameter space for 
the occurrence of several liquid crystalline phases includ- 
ing nematic, columnar, smectic-A, plastic crystalline, and 
orientationally ordered crystalline phases. Recent numer- 
ical studies [47 of the liquid crystalline PFC model in two 
spatial dimensions have shown that a variety of phase co- 
existences occur as a function of the model parameters. 
Therefore, the liquid crystalline PFC model [45 provides 
a simple and direct avenue to access the interface struc- 
ture, which still incorporates the correct physics. 

As a result, we find that the width of the interface 
with respect to the nematic order parameter is larger 
than the width of the mean density interface. In ap- 
proaching the interface from the plastic crystalline side, 
at first the mean density goes down and then the nematic 
order parameter follows. The relative shift of these two 
profiles can be larger than a full lattice constant of the 
plastic crystal. Finally, we also present numerical results 
for the dynamic relaxation of an initial order-parameter 
profile towards its equilibrium inter facial profile. Our 
results can in principle be verified in real-space experi- 
ments of colloidal dispersions, which can be confined to 
monolayers [ 48ri51j . A transient non-monotonic behav- 
ior of the conserved mean-density profiles occurs, which 
is much more pronounced than non-monotonicities in the 
non-conserved orient ational order profile. 

The paper is organized as follows: after the presenta- 
tion of a suitable PFC model for liquid crystals in Sec. 
\n\ we describe a numerical method for the solution of 



the dynamical PFC equations in Sec. IIP Our results ob- 



tained by numerical calculations are discussed in Sec. [IV 
Finally, we conclude in Sec. [V| 



tensor can be parametrized as 

Qij{f,t) = S{r,t)(ni{r,t)nj{r,t) - (1) 

with the nematic order parameter S{f,t) and the (nor- 
malized) nematic director n{f^t) = (ni,n2) (see Refs. 
[Ill [111 [52]). 

A. Static free-energy functional 

The static properties of a liquid crystalline system are 
described by a free-energy functional Qij], which is 
minimized with respect to tp{r) and Qij{r) in thermo- 
dynamic equilibrium. After an appropriate rescaling of 
the length and energy scales, this free-energy functional 
obtains the dimensionless form^ [53] 
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3 6^^ ^4 



(2) 



B3{dkt^){diQki) + DiQli + D2{diQ 



kl) 



with the Laplace operator A = 9| and the five dimen- 
sionless coupling parameters Ai, A2, ^3, I^i, and 1^2 • 



B. Dynamical equations 

The corresponding dynamical equations of '^(r, t) and 
Qij (^, ^) can be derived from classical dynamical density 
functional theory [54 and are given by [53] 



(3) 
(4) 



II. PFC MODEL FOR LIQUID CRYSTALS IN 
TWO SPATIAL DIMENSIONS 



A PFC model for apolar^ liquid crystals in two spatial 
dimensions was given in Refs. [42l |45l [47l [52l [53]. It 

describes the static properties and dynamical behavior of 
a liquid crystalline system in terms of two dimensionless 
order-parameter fields: the reduced translational density 
^/^(r, t) and the symmetric and traceless nematic tensor 
Qij{r,t) with position r = (x^y) and time t. For liquid 
crystalline particles with a symmetry axis, the nematic 



with the dimensionless current jf (r , t) and the dimen- 
sionless quasi-current <l>^(r, t). In constant-mobility ap- 
proximation, this current and quasi-current are given by 



- 

^9 = 



(5) 



(6) 



^ We neglect a possible macroscopic polarization. 



^ Einstein's sum convention is used throughout this paper. No- 
tice that powers of indexed quantities involve repeated indices 
and thus summation, i.e., for example, Q^j = Qij Qij = 

Qij Qij- 
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with the three dimensionless mobihty parameters ai, 0^3, 
and 0^4 and the thermodynamic conjugates 



5F_ 



(5 J- 



(7) 



of ?/^(r, t) and Qij{f^ t), respectively. The thermodynamic 
conjugates fohow directly from the free-energy functional 
([2| by functional differentiation: 

^^=-^'^ + (2V^ - 1)^ + 2A,^ 

- 2D2dk{diQkj + S^Q/ci - ^ij{diQki)) ' 

For a comparison of the dimensionless rescaled param- 
eters in Eqs. ([2|, ([5|, and (|6| with the corresponding 
parameters in the notation of Refs. [42l |47l see ap- 
pendix [A| 



III. NUMERICAL SOLUTION OF THE PFC 
MODEL 

By inserting Eqs. ([8| and (|9| into Eqs. ([5| and (|6| one 
obtains for the dynamics (|3| and Q a system of six cou- 
pled nonlinear partial differential equations. In order to 
solve this system numerically, we decouple and linearize 
it. A simplification is possible due to the symmetry and 
tracelessness of the nematic tensor. Defining the vari- 
ables Qi = Qi^i and = 1, we can write the system 
of dynamical equations as 



(10) 



with the operator ▲ = — 92^2, 28182) that is related 
to the Cauchy-Riemann operator. The thermodynamic 
conjugates reformulated in the new variables read 

with the polynomials 



(12) 



Tn = — we obtain (by treating some terms ex- 
plicitly and others implicitly) the decoupled systems 



+ 1 



2asAiq, 



V^1+i-2(Ai-A2(A + A2))V^,+i 

= CJ^('0n+l, ^n) — Bsk^qi^n 



(13) 
(14) 



4(aiA - 2^4 



— + 2^3 A^V^, 



and 

Qi,n-\-l 

<7i,n+l - 2(2i:>i - D2A)qi^ri+l = UJqi^Pn, Qn+l) 

- BsAilpn . 



(15) 



Linearizing uj^{ilj^q) and ujq{il)^q)i around the old time 
step tn^ two linear systems can be solved one after the 
other for all n. The linearizations of the polynomials ( 12 ) 
read 

^^(V^n, ^n) = (- 2?An + 2?/^^ + q^^) 



UJq{^lJn,qn,qn+l)i 



1 



{Ai^n{lpn-^) + qln) 



(16) 



1 



V Qi,nQj ,nQj ,n-\-l 



However, such a simple time stepping scheme with con- 
stant step size turns out to be impractical, if stationary 
configurations have to be obtained in the simulations. 
We therefore used a higher-order embedded Rosenbrock 
scheme with an adequate step size control for the time 
discretization. A detailed description of this scheme con- 
cerning some numerical issues will be given elsewhere. 
For the discretization in space we used the finite element 
method. 



IV. RESULTS 

We first restrict ourselves to certain parameter com- 
binations, which allow for several liquid crystalline co- 
existences. In detail, we fix the parameters A2 = 14, 
Bs — —0.4, = 1, and D2 = 0.8, but vary the param- 
eter Ai (which corresponds to some formal temperature 
in the context of mean-field theories) and the reduced 
mean density tp.^ The resulting equilibrium bulk phase 
diagram is shown in Fig. [T] in consistency with earlier 
data [4T. In the parameter range of Ai and tjj shown, 
the phase diagram exhibits three stable liquid crystalline 
phases, namely, the isotropic phase, a plastic triangular 



We could discretize in time using a semi-implicit Euler 
discretization. Let ^1,^2,^3... be a sequence of time 
steps. Defining V^(^,U), Qi^n qiif.tn), and 



The parameters in the dynamical equations (|3|-(|6| are always 
chosen to be q;i = 0:3 = 0:4 = 1. Clearly, the stationary results 
do not depend on their particular values. 
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FIG. 1: Phase diagram with coexistence regions for the mean 
density -0 G [—1.6,0.4] and the parameters Ai G [1.3,3.5], 
A2 14, B3 -0.4, Di 1, and D2 0.8. Three 
different hquid crystalhne phases are reahzed: isotropic, 
columnar/smectic A (CSA), and plastic triangular crystalline 
(PTC). The coexistence regions (shaded areas) are calcu- 
lated using Maxwell's double tangent construction. The black 
dashed lines in the coexistence regions indicate the intersec- 
tion lines of the energy curves of the two adjacent phases. 
Six black circles indicate certain parameter combinations for 
which detailed calculations were performed (see Figs, [s] [7| 
and [4][6t . 



crystal (PTC)^, and a columnar phase. As we consider 
two spatial dimensions here, a columnar phase is indis- 
tinguishable from a smectic A phase, therefore we call 
the latter columnar/smectic A (CSA)^ phase. The coex- 
istence regions, as obtained by a Maxwell double tangent 
construction, are depicted by the shaded area in Fig. [T] 
We selected in total six different coexistence conditions 
as labeled by black circles in Fig. [l] which correspond 
to three isotropic-PTC and three CSA-PTC coexistence 
situations serving as basic reference situations for our 
subsequent investigations. 

A typical example for an isotropic-PTC interfacial pro- 
file is presented in Fig. |2]for the (lO)-orientation of the 
hexagonal crystal. In the bulk PTC phase, there are peri- 
odic peaks in the full density profile ip{x,y) at the crystal 
lattice positions, shown as a contour plot in Fig. [2j The 
typical standard deviation of these peaks (the so-called 
Lindemann parameter) is pretty large with about 27% of 
the lattice constant. The corresponding orient ational or- 
dering as embodied in the nematic tensor is complicated 
and exhibits topological defects in the Wigner-Seitz cell 



The plastic triangular crystal in phase diagram[l]is called "plastic 
triangular crystal 2" (PTC2) in Ref. 4T. 
^ This columnar/smectic A (CSA) phase is called "C/SA phase" 
in Ref. l47l. 
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FIG. 2: Top: two contour plots for ip(x,y) and S{x,y) at 
an isotropic-PTC coexistence with Ai = 3.21 and = —0.3 
(the other parameters are the same as in Fig.[T]). n{r) is rep- 
resented by short black lines that are superimposed to the 
lower contour plot. Bottom: averaged density {ip){x) (left or- 
dinate) and averaged nematic order parameter (S){x) (right 
ordinate). The x-direction is chosen perpendicular to the in- 
terface, while the ^-axis is parallel to the interface. The aver- 
aged quantities are defined by (/)(x) = f dy' f^^^dx' f{x\ y') 

for / G {V^, S} with the width of the stripes 2a — 47r/(/cv^) 
and k = l/\/2. 



of the lattice, see Refs. [47l |55] for a more detailed dis- 
cussion. The mean orientational unit vector field n{x^y) 
as obtained by the direction of the eigenvector of the 
nematic tensor corresponding to the highest eigenvalue, 
is sketched by short black lines in Fig. |2j The largest 
eigenvalue itself, the scalar nematic order parameter field 
S{x^y), is also presented as a contour plot in Fig. [2] In 
the isotropic phase, on the other hand, the density field 
is constant and the nematic order parameter vanishes. 
In between there is an interfacial region with laterally 
averaged profiles {il^){x) and {S){x) with x denoting the 
direction perpendicular to the interface (see caption of 
Fig.[2|. 

We define a typical interface width of an order pa- 
rameter profile f{x^y) G {V^(x, ^), 6'(x, ?/)} as the dis- 
tance of the positions, where a tanh- approximation of 
{f){x) attains the values 0.95(/)(-oo) + 0.05(/)(oo) and 
0.05(/)(— 00) + 0.95(/)(oo), respectively. These widths 
for ijj^x^y) and S{x^y) are indicated in Fig. [2] Re- 
markably, the width of the density profile is signifi- 
cantly smaller than the width of the orientational pro- 
file. The position, where the tanh- approximation of an 
averaged field {f){x) with / G {i^^S} attains the value 
((/)(— 00) + (/)(oo))/2 can be taken as a natural location 
^(/) of the interface with respect to this field. Interest- 
ingly, as revealed in Fig. |2j the location of the averaged 
density profile {ip){x) and the averaged orientational pro- 
file {S){x) do not coincide. The location of the orienta- 
tional profile is more shifted towards the isotropic phase 
than the location of the density profile. This means that 



coming from the isotropic side, at first the nematic or- 
der builds up and then the density fohows. This finding 
is reminiscent to the fiuid-crystal interface of systems of 
spherical particles [23, 56 , which can be described by a 
two-order-parameter description involving the conserved 
mean density and a non-conserved crystallinity [TTl [57] . 
Coming from the fluid side, also in this case, the non- 
conserved crystallinity starts to grow first and the density 
follows. 

We have further studied the dependence of the inter- 
face widths on the parameter Ai. As Ai is increased, 
the coexistence comes closer to a critical point where the 
interfacial widths diverge. This trend is documented in 
Fig. [3] The upper plot in Fig. [3] also shows that the width 

Isotropic-PTC interface 
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FIG. 3: Widths of the interfaces of ip{r) and S{r) 
in dependence of Ai. The parameters are {ip,Ai) G 
{(-0.05, 1.6), (0.05, 2.4), (0.31, 3.5)} and for the rest as in Fig. 
[l] Notice that the presented data correspond to the six points 
highlighted by black circles in Fig. [l] and that they are con- 
nected by spline interpolation. In the lower plot the stripes 
of the CSA phase are oriented perpendicular to the interface 
(see Fig. [5]). 

of the orientational order-parameter profile is larger than 
that of the density interface over the full range of Ai . All 
trends are the same for different parameter combinations 
for the isotropic-PTC interface, as documented by Fig. 

El 

Next, we consider the coexistence between the PTC 
and the CSA phase. In this case, the interface structure 
depends on the relative orientations of the two phases. 
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Isotropic-PTC interface 
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FIG. 4: Interface of the isotropic-PTC phase coexistence for 
the same parameters as in Fig. [l] The plots show the trans- 
lational density ip(r) and the nematic order parameter S(f) 
both for a large area and for a close-up view of the interface, 
where blue and red indicate low and high values, respectively. 
In addition, the director field n(r) is represented by short 
black lines that are superimposed to the lowest plots. 



While we fix the orientation of the PTC phase in the 
(lO)-direction, we consider here two possibilities of the 
column direction relative to the interface, namely per- 
pendicular and parallel. For these two different relative 
orientations, the order-parameter fields are given in Figs. 
|5] and [6] for two different parameter combinations of co- 
existence. For perpendicular column direction (see Fig. 



CSA-PTC interface (stripes perpendicular to interface) 











tjj = -0.05, Ai = 1.6 


tp = 0.31, Ai = 3.5 



FIG. 5: The same as in Fig. [i] but now for the CSA-PTC 
interface. Note that the stripes of the CSA phase are oriented 
perpendicular to the interface. 



[5|, the density field reveals that the columns end at a 
lattice density peak. This implies that the degeneracy 
of the column positions is broken by the presence of the 
crystal, which pins the transversal columnar order by the 
interface. Along the columns away from the interface, 
there are still some density undulations in x-direction. 
For parallel column direction (see Fig. [6|, on the other 
hand, there is a nontrivial density field across the inter- 
face insofar as the columns are significantly bent in the 
presence of the crystalline peaks, i. e., the crystal induces 
a systematic undulation of the neighboring columns. The 
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CSA-PTC interface (stripes parallel to interface) 




FIG. 6: The same as in Fig. |5j but now for a CSA-PTC 
interface, where the stripes of the CSA phase are oriented 
parallel to the interface. 

amplitude of this undulation decreases farer away from 
the interface position. Likewise, along the columns there 
is a periodic density modulation in ^/-direction induced 
by the crystalline peaks nearby. 

Results for the interfacial width, similarly defined as in 
the previous case, are shown in the lower plot in Fig. |3j 
where the same trends are observed as for the isotropic- 
PTC interface [see the upper plot in Fig. [s]. The width 
of the orientational interface is considerably larger than 
that for the density profile and there is a strong depen- 
dence on the parameter Ai with huge interfacial widths, 
where the parameters are close to criticality. Like the 
isotropic-PTC interface, the interface position of the den- 
sity profile is more in the PTC-phase than the interface 
position of the orientational profile, which is more in the 
coexisting CSA phase (see Fig. ItI). The shift in the two 




I — ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 1 

1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 

FIG. 7: Distance A^(^, S) =^(3)- ^(^) of the interfaces of 
ip{f) and S{f) in dependence of (we always consider the 
transition from the crystalline to the non-crystalline phase). 
The parameters are the same as in Fig.|3]and the stripes of the 
CSA phase are again oriented perpendicular to the interface 
(see Fig. [5]). 

interface positions depends on the parameters, as shown 
in Fig. [7l 

Finally, we show some results on the dynamical evo- 
lution of the interfacial profiles based on the physical 



dynamics described in Eqs. (|3| and Q. It is important 
to note that the density is a conserved order parameter, 
while the nematic ordering is non-conserved. We plot an 
example of the interface relaxation towards equilibrium 
for a prescribed starting profile in Fig. [8) The orienta- 
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FIG. 8: Time evolution of the averaged order parameters 
('0)(x) (top) and {S){x) (bottom) for an isotropic-PTC co- 
existence. The parameters are ai = as — cx^ — 1, —0.3, 
A\ — 3.21 and for the rest as in Fig. ^ Snapshots are taken 
at times t = 100, 500, 3000, and 10000. At time t = the av- 
eraged translational density is constant {{ip){x) = —0.3) and 
the averaged nematic order parameter {S){x) is a smeared 
Heaviside step function. 



tional order-parameter field is a smeared Heaviside step 
function, while the density is constant. Similar set-ups 
for interfacial kinetics have been studied earlier [58 . The 
density field subsequently takes up the orientational inho- 
mogeneity and both order parameters relax to their equi- 
librium profiles. The density develops a marked transient 
non-monotonic profile and relaxes much slower than the 
orientational order. It takes quite a long time in units of 
the basic time scale of the dimensionless dynamical equa- 
tions ([3|-([6| to end up in the final equilibrium state. Our 
finding shows that in principle our dynamical equations 
(|3| and Q, which reflect the diffusive dynamics of col- 
loidal systems, can be applied to plenty of further growth 
phenomena in the future, which are, however, beyond the 
scope of the present paper. 
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V. CONCLUSIONS 

In conclusion, we have explored the equilibrium struc- 
ture of interfaces between various coexisting liquid crys- 
talline phases using a PFC model for liquid crystals. 
In two spatial dimensions, we have considered explic- 
itly the isotropic-plastic crystalline and the smectic A- 
plastic crystalline interface, which are both anisotropic, 
i.e., they depend on the relative orientation of the two 
coexisting phases. To determine the equilibrium struc- 
tures numerically, we calculated the relaxation of the 
dissipative PFC dynamics towards equilibrium (i.e., the 
minimization of the PFC functional) under the constant- 
mobility approximation using the finite element method. 

Basically, we have considered a two-order-parameter 
description of the interfaces containing the conserved 
(translational) density field and the non-conserved (ori- 
entational) nematic tensor. The phase diagram, the typi- 
cal widths of the interfaces, the order-parameter profiles, 
and their dynamics were computed. For the isotropic- 
plastic crystalline interface we find that in approaching 
the interface from the isotropic side, at first the nematic 
order builds up and then the density follows. The rel- 
ative shift of the two profiles is about half the lattice 
constant of the plastic crystal. This finding is reminis- 
cent to the fiuid-crystal interface of systems of spherical 
particles ^23i i56| , which can be described by a two-order- 
parameter description involving the conserved mean den- 
sity and a non-conserved crystallinity [HI [57j . For the 
fiuid-crystal interface, a similar shift has been found: if 
the interface is approached from the fiuid side, first the 
(non-conserved) crystallinity increases and then the (con- 
served) mean density follows [JT, .59-61^ . This has to do 
with the fact that a fiuid is more responsive to an os- 
cillatory density wave than to a global density change 

Our results can be verified either in particle-resolved 
computer simulations [62^ or in experiments. Particle- 
resolved computer simulations for rod-like systems have 
been performed both for structure [4l I63H65] and dy- 
namics j66l [67] in various situations. So far experi- 
ments are concerned, most notably colloidal liquid crys- 
tals [niEHlES] that are confined to two spatial dimensions 
are ideal realizations of our model. One important ex- 
ample is a suspension of the tobacco mosaic virus, which 
can be confined to monolayers and which shows a va- 
riety of liquid crystalline phases [8] , but there are more 
other examples of liquid crystalline rod-like particle sus- 
pensions, which have been prepared in a controlled way 
(see, e.g., Refs. [48, 49, 51 ). 

Future work should extend the present study to three 
spatial dimensions |42l |46], which would require more 
numerical work but promises a richer equilibrium bulk 
phase diagram. Also the dynamics of a growing crys- 
talline front, which has been studied for spherical par- 
ticles already in detail [7QH72] . should be addressed for 
liquid crystals as well. If a plastic crystalline phase grows 
into an isotropic phase, it would be interesting to fol- 



low the origin of topological defects in the director field, 
which have to grow out of nothing. Moreover, crystal- 
fiuid interfaces in external fields, like gravity, exhibit un- 
usual effects already for isotropic particles [73"^ and it 
would be challenging to explore this for liquid crystalline 
interfaces [75 . Finally, our model should be generalized 
to liquid crystals on manifolds [76 to describe nematic 
[77] or smectic bubbles ^78j . 



Appendix A: Notation 

Since the PFC model presented in Sec. [IT] is equivalent 
to PFC models given in different notation in Refs. [42l 
1471 f53] , we here clarify the relationship of our notation 
to the notation used in the literature. This especially 
simplifies the comparison of our phase diagram [l] to the 
corresponding phase diagrams in Ref. [47 . 

If we denote the eight parameters in Eqs. d2|, ([5|, and 
(|6| with a prime (i.e., ^3, a^, 0^3, 0^4) to 

avoid confusion with a similar notation in Refs. [42] [53], 
the characteristic length and the characteristic energy 
E'^^ which have been chosen to make the PFC model in 
the pr esent arti cle dimensionless, can be expressed by 
Ic = \J-MIM and E'^ = -{irp/ (3){As/A2) in terms of 
the reference particle number density p, the inverse ther- 
mal energy /3, and the parameters A2 and As in Ref. 
[53] , In Ref. [42], the notation is analogously, but with 
Pref instead of p. 

The parameters in the free-energy functional ([2| can 
be related to the parameters in Refs. [42] [53j by 



and 



27rp ' 



A' 



^2 



27rpA3 



B' 



AoDo 



A2BS 
27rpAs 



(Al) 



(A2) 



4 27rp ' ^ 27TpAs 

In case of the current jsj and the quasi-current (|6|, a 
comparison with Refs. |42l [53] leads to the relations 



tr.Ec 



a. 



ao 



IT 



as , a'^ = tcE^a^ . (A3) 



In Ref. [47 , a different notation is used. A comparison 
of the free-energy functional Q with the corresponding 
free-energy functional in Ref. |47 leads to 



A[ = Bi , = 45. , B's 



-4F 



(A4) 



and 



D[ = 2D , D'2 = SE . 



(A5) 

Furthermore, the length and energy scales of these two 
free-energy functionals are different. If and E^^ denote 
the characteristic length and energy in the present article 
and Ic and Ec denote the corresponding quantities in Ref. 
[47 , they can be related to each other by 



K = Ie, 

c 2 



(A6) 
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